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Abstract 

We will display the fundamental structure of classical electrody- 
namics. Starting from the axioms of (1) electric charge conservation, 
(2) the existence of a Lorentz force density, and (3) magnetic flux 
conservation, we will derive Maxwell's equations. They are expressed 
in terms of the field strengths (E,B), the excitations (T>,H), and the 
sources (p,j). This fundamental set of four microphysical equations 
has to be supplemented by somewhat less general constitutive assump- 
tions in order to make it a fully determined system with a well-posed 
initial value problem. It is only at this stage that a distance concept 
(metric) is required for space-time. We will discuss one set of possible 
constitutive assumptions, namely T> ~ E and H ~ B. file erik8a.tex, 
1999-07-27 

* Dedication to Erik W. Grafarend on the occasion of his 60th birthday: Wir wissen, 
dass heutzutage auch die Geodaten relativistische Effekte bei ihren "Triangulationen" 
beriicksichtigen miissen. Wohl auch aus diesem Grunde hat Herr Grafarend immer ein 
ofFcnes Ohr fur entsprechende Theorien gehabt. Die einfachste relativistische klassische 
Feldtheorie, die wir kennen, ist die Elektrodynamik. Wir widmen diese Ausarbeitung 
Herrn Grafarend zu seinem 60.Geburtstage in der Hoffhung, dass er sich iiber die schonen 
Seiten dieser Darstellung genauso freut, wie wir es tun. Und dies umso mehr, als dass 
Herr Grafarend gleich am Anfang seiner Karriere sich intensiv mit Geometrie und Cartan- 
Formalismus auseinandergesetzt hat. 

^Permanent address: Moscow State University. 
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1 Introduction 



Is it worthwhile to reinvent classical electrodynamics after it has been with 
us for more than a century? And after its quantized version, quantum elec- 
trodynamics (unified with the weak interaction) had turned out to be one of 
the most accurately tested successful theories? We believe that the answer 
should be affirmative. Moreover, we believe that this reformulation should 
be done such that it is also comprehensible and useful for experimental physi- 
cists and (electrical) engineers^. 

Let us collect some of the reasons in favor of such a reformulation. First 
of all an "axiomatics" of electrodynamics should allow us to make the funda- 



mental structure of electrodynamics transparent, see, e.g., Sommerfeld (16 
or (IJ, 0, |2"0| . We will follow the tradition of Kottler-Cartan-van Dantzig, see 
Truesdell & Toupin and Post [|10j, and base our theory on two experi- 
mentally well established axioms expressed in terms of integrals, conservation 
of electric charge and magnetic flux, and a local axiom, the existence of the 
Lorentz force. All three axioms can be formulated in a 4-dimensional (space- 
time) continuum without using the distance concept (i.e. without the use of a 
metric), see Schrodinger |I5| |. Only the fourth axiom, a suitable constitutive 



law, is specific for the "material" under consideration which is interacting 
with the electromagnetic field. The vacuum is a particular example of such 
a material. In the fourth axiom, the distance concept eventually shows up 
and gives the 4-dimensional continuum an additional structure. 

Some of the questions one can answer with the help of such a general 
framework are: Is the electric excitation D a microscopic quantity like the 
field strength El Is it justified to give T> another dimension than El The 
analogous questions can be posed for the magnetic excitation H and the field 
strength B. Should we expect a magnetic monopole and an explicit magnetic 
charge to arise in such an electrodynamic framework? Can we immediately 
pinpoint the (metric-independent) constitutive law for a 2-dimensional elec- 
tron gas in the theory of the quantum Hall effect? Does the non-linear Born- 
Infeld electrodynamics fit into this general scheme? How do Maxwell's equa- 
tions look in an accelerated reference frame or in a strong gravitational field 
as around a neutron star? How do they look in a possible non-Riemannian 
spacetime? Is a possible pseudoscalar axion field compatible with electrody- 



^or this reason, we apply in this article the more widespread formalism of tensor 
analysis ("Ricci calculus", see Schouten [|l4)) rather than that of exterior differential forms 
("Cartan calculus", see Frankel ||) which we basically prefer. 
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namics? And eventually, on a more formal level, is the calculus of exterior 
differential forms more appropriate for describing electrodynamics than the 
3-dimensional Euclidean vector calculus and its 4-dimensional generalization? 
Can the metric of spacetime be derived from suitable assumptions about the 
constitutive law? 

It is really the status of electrodynamics within the whole of physics which 
comes much clearer into focus if one follows up such an axiomatic approach. 

2 Foliation of the 4-dimensional spacetime con- 
tinuum 

From a modern relativistic point of view, the formulation of electrodynam- 
ics has to take place in a 4-dimensional continuum (differentiable manifold) 
which eventually is to be identified with spacetime, i.e. with a continuum de- 
scribed by one "time" coordinate x° and three "space" coordinates X - X ■ X 
or, in short, by coordinates x\ with i — 0, 1, 2, 3. Let us suppress one space 
dimension in order to be able to depict the 4-dimensional as a 3-dimensional 
continuum, as shown in Fig.Q. 




Figure 1: Foliation of spacetime: Each hypersurface h a represents, at a time 
a, the 3-dimensional space of our perception one dimension of which is sup- 
pressed in the figure. The positive time direction runs upwards. 
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We assume that this continuum admits a foliation into a succession of dif- 
ferent leaves or hypersurfaces h. Accordingly, spacetime looks like a pile of 
leaves which can be numbered by a monotonically increasing (time) parame- 
ter a. A leaf h a is defined by a{x l ) = const. It represents, at a certain time a, 
the ordinary 3-dimensional space surrounding us (in Fig.|l| it is 2-dimensional, 
since one dimension is suppressed). 

At any given point in h a , we can introduce the covector fcj := <9jcx and 
a 4-vector n = (ri 1 ) = (n , n 1 , n 2 , n 3 ) = (n°,n a ) such that n is normalized 
according to 

n* h = n%a = 1 . (1) 

Here a,b . . . = 1, 2, 3 and = 0,1, 2, 3. Furthermore, summation over 

repeated indices is always understood. The vector n l is "normal" to the leaf 
h a , whereas the covector ki is tangential^ to it. 

With the pair (n, k) we can construct projectors which decompose all 
tensor quantities into longitudinal and transversal constituents with respect 
to the vector n, see Fig.||. Indeed, the matrices 

Vj := n { kj and Tj := 8) - n* fej with V$ + Tj = 5) , (2) 

represent projection operators, i.e. 

Vj V k = L\ , T l j T 3 k = T\ , Vj T 3 k = T l j V k = . (3) 

Taking an arbitrary covector Ui, we now can write it as 

Ui = L Ui + Ui , where ± U l := L\ Uj and := T j t Uj . (4) 

Obviously Ui describes the longitudinal component of the covector and XV 
its transversal component, with n % XV = 0. Analogously, for an arbitrary 
vector V 1 , we can write 

V 1 = V* + Yl , where ^ := Vj V j and F := V j V j . (5) 

Its transversal component V_ l fulfills V? ki = 0. This pattern can be straight- 
forwardly generalized to all tensorial quantities of spacetime. 

For simplicity, we confine our attention to the particular case when "adapt- 
ed" coordinates x % = (a, x a ) are used and when the "spatial" components of 
n vanish, i.e., n % = (1,0,0,0). In that case, we simply have ki = diO = 
(1, 0, 0, 0) and hence a can be treated as a formal "time" coordinate. 

2 The term "tangential" is used here in the sense of exterior calculus in which a covector 
(or 1-form) is represented by two ordered parallel planes - and the first plane is tangential 
to h a . 
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3 Conservation of electric charge (axiom 1) 



The conservation of electric charge was already recognized as fundamental 
law during the time of Franklin (around 1750) well before Coulomb discov- 
ered the force law in 1785. Nowadays, at a time, at which one can catch 
single electrons and single protons in traps and can count them individu- 
ally, we are more sure than ever that electric charge conservation is a valid 
fundamental law of nature. Therefore matter carries as a primary quality 
something called electric charge which only occurs in positive or negative 
units of an elementary charge e (or, in the case of quarks, in l/3th of it) 
and which can be counted in principle. Thus it is justified to introduce the 
physical dimension of charge q as a new and independent concept. Ideally 
one should measure a charge in units of e/3. However, for practical reasons, 
the Si-unit C (Coulomb) is used in laboratory physics. 




longitudinal to n and the transversal piece J_, global conservation of charge. 

Two remarks are in order: Charge is an additive (or extensive) quantity 
that characterizes the source of the electromagnetic field. It is prior to the 
notion of the electric field strength. Therefore it is not reasonable to measure, 
as is done in the CGS-system of units, the additive quantity charge in terms 
of the unit of force by applying Coulomb's law. Coulomb's law has no direct 
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relation to charge conservation. Secondly, in the Si-system, for reasons of 
better realization, the Ampere A as current is chosen as the new fundamental 
unit rather than the Coulomb. We have C = As (s = second). 

As a preliminary step, let us remind ourselves that, in a 4-dimensional 
picture, the motion of a point particle is described, as in Fig.|], by a curve 
in spacetime, by a so-called worldline. The tangent vectors of this worldline 
represent the 4-velocity of the particle. 

If we mark a 3-dimensional volume which belong to a certain hyper- 
surface h/j , then the total electric charge inside Q3 is 



with p as the electric charge density. The total charge in space, which we find 
by integration over the whole of space, i.e., by letting Q 3 — > h a , is globally 
conserved. Therefore the integral in (|^) over each hypersurface h ai , h a2 , . . . 
keeps the same value. 

The local conservation of charge, see Fig.|3|, translates into the following 
fact: If a number of worldlines of particles with one elementary charge enter a 
prescribed but arbitrary 4-dimensional volume f2 4 , then, in classical physics, 
the same number has to leave the volume. If we count the entering worldlines 
as negative and the leaving ones as positive (in conformity with the direction 
of their normal vectors), then the (3-dimensional) surface integral over the 
number of worldlines has to vanish. 

Now, the natural extensive quantities to be integrated over a 3-dimen- 
sional hypersurface are vector densities, see the appendix. Accordingly, in 
nature there should exist a 4-vector density J % with 4 independent compo- 
nents which measures the charge piercing through an arbitrary 3-dimensional 
hypersurface. Therefore, it generalizes in a consistent 4-dimensional formal- 
ism the familiar concepts of charge density p and current density j a . The 
axiom of local charge conservation then reads 



where the integral is taken over the (3-dimensional) boundary of an arbitrary 
4-dimensional volume of spacetime, with d 3 Si being the 3-surface element, 
as defined in the appendix. 




(6) 




(7) 
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Figure 3: Local conservation of charge: Each worldline of a charged particle 
that enters the finite 4- volume f2 4 via its boundary dfl^ has also to leave f^. 

If we apply Stokes' theorem, then we can transform the 3-surface integral 
in (0) into a 4-volume integral: 

J(d t T)d 4 S = 0. (8) 

Since this is valid for an arbitrary 4-volume f2 4 , we find the local version of 
the charge conservation as 

diT = . (9) 

In this form, the law of conservation of charge is valid in arbitrary coordi- 
nates. 

If one defines a particular foliation, then one can rewrite (||) in terms 
of decomposed quantities that are longitudinal and transversal to the corre- 
sponding normal vector n. The 4-vector density J % decomposes as 

ji = ±ji + ji _ ( 1Q ) 

When adapted coordinates are used, the decomposition procedure simplifies 
and allows to define the 3-dimensional densities of charge p and of current 

j a as 
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With this, one can rewrite the definition of charge (|G|) in an explicitly coor- 
dinate invariant form 

Q= J J* d 3 Si , (12) 

since on h a we have d 3 So = dx 1 dx 2 dx 3 and d 3 S a = 0. Furthermore, Eq.(|9]) 
can be rewritten in (l+3)-form as the more familiar continuity equation 

d a p + d a j a = 0. (13) 

The charge Q in fll2"|) has the absolute dimension^ q. The 4-current is a 
density in spacetime, and we have [J'} = q/{tl 3 ). Thus the components carry 
the dimensions [p] = [J°] = q/l 3 = C/m 3 and [j a ] = [J a ] = q/{tl 2 ) = A/m 2 . 



4 The inhomogeneous Maxwell equations as 
consequence 

Because of axiom 1 and according to a theorem of de Rham, see ||, the 
electric current density from (0) or (H) can be represented as a "divergence" 
of the electromagnetic excitation: 

X = dj n ij , K j = -n ji . (14) 

The excitation Ti^ is a contravariant antisymmetric tensor density and has 
6 independent components. One can verify that, due to the antisymmetry of 
7i %: > , the conservation law is automatically fulfilled, i.e., diJ 1 = didfH,^ = 0. 

The 4-dimensional set ([14]) represents the inhomogeneous Maxwell equa- 
tions. They surface here in a very natural way as a result of charge con- 
servation. Charge conservation should not be looked at as a consequence of 
the inhomogeneous Maxwell equations, but rather the other way round, as 
shown in this tutorial. Of course, 7i JJ is not yet fully determined since 

m =H ij + e ijkl d k ^ L (15) 

3 A theory of dimensions, which we are using, can be found in Post [[To) , e.g.. A quantity 
has an absolute dimension, and if it is a density in spacetime we divide by tl 3 . The 
components pick up a t (a t" 1 ) for an upper (a lower) temporal index and an I (an Z -1 ) for 
an upper (a lower) spatial index. A statement, see |4||, that E and B must have the same 
dimension since they transform into each other is empty without specifying the underlying 
theory of dimensions. 
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also satisfies ( fR|) for an arbitrary covector field ipi- 

The (1 + 3)-decomposition of H,^ is obtained similarly to the decomposi- 
tion of the current (|I0|): 

H ij = ± H ij +H ij . (16) 
The nontrivial components of the longitudinal and transversal parts read 

H 0a = ± H 0a = pa^ ^afe = yob = fbc R ^ ^ 

with the electric excitation T> a (historical name: "dielectric displacement") 
and the magnetic excitation H a ("magnetic field"). Here e abc is the totally 
antisymmetric 3-dimensional Levi-Civita tensor density with e 123 = 1. 

If we substitute the decompositions (|10D and (|16D into QUI) , we recover 



the 3-dimensional form of the inhomogeneous Maxwell equations, 

d a V a = p, e abc d b H c -d a V a =j\ (18) 

or, in symbolic notation, 

div£> = p, cm\H-V = j. (19) 

Since electric charge conservation is valid in microphysics, the corresponding 
Maxwell equations ( |T8"D or flT5| ) are also microphysical equations and with 
them the excitations T> a and H a are microphysical quantities likewise - in 
contrast to what is stated in most textbooks, see and \T2\, compare also 
0,e.g.. 

From ( [IB] ) we can immediately read off [D a ] = [I p] = q/l 2 = Cjvn? 

SI 

and [H a ] = [If] = q/(tl) = A/m. Before we ever talked about forces 
on charges, charge conservation alone gave us the inhomogeneous Maxwell 
equations including the appropriate dimensions for the excitations T> a and 
H a . 

Under the assumption that T> a vanishes inside an ideal electric conductor, 
one can get rid of the indeterminacy of T> a , as spelled out in flTo]), and we can 
measure T> a by means of two identical conducting plates ("Maxwellian dou- 
ble plates" ) which touch each other and which are separated in the T> a - field 
to be measured. The charge on one plate is then measured. Analogous re- 
marks apply to H a . Accordingly, the excitations do have a direct operational 
significance. 
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5 Force and field strengths (axiom 2) 



By now we have exhausted the information contained in the axiom 1 of 
charge conservation. We have to introduce new concepts in order to com- 
plete the fundamental structure of Maxwell's theory. Whereas the excita- 
tion Ti. = (V a , H a ) is linked to the charge current J = (p, j a ), the electric 
and magnetic field strengths are usually introduced as forces acting on unit 
charges at rest or in motion, respectively. In the purely electric case with a 
test charge q, we have in terms of components 

F a ~qE a , (20) 

with F as force and E as electric field covector. 

Let us take a (delta-function-like) test charge current J = (p, j a ) centered 
around a point with spatial coordinates x a . Generalizing (|20|) , the simplest 
relativistic ansatz for defining the electromagnetic field reads: 

force density ~ field strength x charge current density . (21) 

We know from Lagrangian mechanics that the force ~ dL/dx l is represented 
by a covector with the absolute dimension of action h (here h is not the 
Planck constant but rather only denotes its dimension). Accordingly, with 
the covectorial force density f i: the ansatz ( pTj) can be made more precise as 
axiom 2: 

fi = Fij Cf^ i Fij = Fji . (22) 

The newly introduced covariant 2nd-rank 4-tensor is the electromagnetic 
field strength. The force density fi was postulated to be normal to the 
current, fiJ % = 0. Thus the antisymmetry of the electromagnetic field 
strength is found, i.e., F^ depends on 6 independent components. We know 
the notion of force from mechanics, the current density we know from axiom 
1. Accordingly, axiom 2 is to be understood as an operational definition of 
the electromagnetic field strength F^. 
With the decomposition 

F« = - L Fy + £ y , (23) 

we find the identifications for the electric field strength E a and the magnetic 
field strength B a (historical names: "magnetic induction" or "magnetic flux 
density" ) : 

F a o = F a0 = E a , F ab = F_ ab = e a bc B c . (24) 
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These identifications are reasonable since for the spatial components of (22) 
we recover the Lorentz force density and, for the static case, Eg. (f2~0|) : 

f a = F aj Ji = F a0 J° + F ab J b = pE a + e abc j b B c . (25) 

Symbolically, we have 

f = pE+jxB. (26) 
The time component of (|22| ) represents the electromagnetic power density: 

fo = F 0a J a = -E a f. (27) 

6 Conservation of magnetic flux (axiom 3) 

Axiom 2 on the Lorentz force gave us a new quantity, the electromagnetic 
field strength with the dimension [F] = action/ charge = h/q =: <f), with 

SI 

<p = work x time /charge = voltage x time = Vs = Wb. Here Wb is the 
abbreviation for Weber. Thus its components carry the following dimensions: 
[E a \ = [F a0 ] = <t>/{tl) = V/m and [B c ] = [F ab ] = 0/Z 2 i Wb/m 2 = T (for 
Tesla) . 

We are in need of an experimentally established law that relates to F. And 
we would prefer, as in the case of the electric charge, to recur to a counting 
procedure. What else can we count in relation to the electromagnetic field? 
Certainly magnetic flux lines in the interior of a type II superconductor which 
is exposed to a sufficiently strong magnetic field. And these flux lines are 
quantized. In fact, they can order in a 2-dimensional triangular Abrikosov 
lattice, see Fig^. These flux lines carry a unit of magnetic flux, a so-called 
flux quantum or fluxon with $ = h/(2e) = 2.07 x 10~ 15 Wb, see Tinkham 
r7||; here h is the Planck constant and e the elementary charge. These flux 



lines can move, via its surface, in or out of the superconductor, but they 
cannot vanish (unless two lines with different sign collide) or spontaneously 
come into existence. In other words, there is a strong experimental evidence 
that magnetic flux is a conserved quantity. 

The number 2 in the relation $ — h/(2e) is due to the fact that the 
Cooper pairs in a superconductor consist of 2 electrons. Moreover, outside a 
superconductor the magnetic flux is not quantized, i.e., we cannot count the 
flux lines there with the same ease that we could use inside. Nevertheless, 
as we shall see, experiments clearly show that the magnetic flux is conserved 
also there. 
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Figure 4: Sketch of an Abrikosov lattice in a type II superconductor in 3-di- 
mensional space. In contrast to all the other figures, this is not an image of 
4-dimensional spacetime. 

As we can take from Fig.[|, the magnetic flux should be defined as a 
2-dimensional spatial integral. These flux lines are additive and we have 

$= J B a d 2 S a . (28) 

Here B a is the magnetic field strength and d 2 S a the spatial 2-surface 
element. This definition of the magnetic flux should be compared with the 
definition (|6|) of the charge. Here, in (|28"D, we integrate only over 2 dimensions 
rather than over 3 dimensions, as in the case of the charge in (|6]). Thus in 
a spacetime picture in which one space dimension is suppressed, see Fig.|5], 
our magnetic flux integral looks like an integral over a finite interval [A, B] 
embedded into the hypersurface h ai . 

Now we are going to argue again as in Sec. 3. If f2 2 — > oo, i.e., if we 
integrate over an infinite spatial 2-surface (A — > +oo, B — > — oo), then the 
total magnetic flux at time <j\ is given by (^). If we propagate that interval 
into the (coordinate) future, see the interval on the hypersurface h a2 , then 
magnetic flux conservation requires the constancy of the integral $. In other 
words, if we orient the integration domain suitably, the loop integral, the 
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Figure 5: Magnetic flux in spacetime: The magnetic field £> a , if integrated 
over the interval [A, B], represents, at time cj, the magnetic flux piercing 
through this 2-dimensional integration domain. 

domain of which is drawn in Fig.|5], has to vanish since no flux is supposed to 
leak out along the dotted "vertical" domains at spatial infinity. 

Analogously as we did in the case of charge conservation, we want to 
formulate a corresponding local conservation law in an explicitly covariant 
way. We saw that the global conservation of magnetic flux is expressed as 
the vanishing of the integral of B over the particular 2-dimensional loop in 
Fig.[5]. In a 4-dimensional covariant formalism, the natural intensive objects 
to be integrated over a 2-dimensional region are second order antisymmetric 
covariant tensors, see the appendix. The magnetic field strength B is just 
a piece of the electromagnetic field strength F. Thus, it is clear that the 
natural local generalization of the magnetic flux conservation, our axiom 3, 
is 

/ ~ F tJ d 2 S^ = , (29) 

where the integral is taken over the boundary of an arbitrary 3-dimensional 
hypersurface of spacetime, as is sketched in Fig|| We apply Stokes' theorem 

J e^dyF^d'S^O, (30) 
n 3 
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Figure 6: Conservation of magnetic flux in spacetime: Consider the arbi- 
trary 3-dimensional integration domain Q 3 . The integral vanishes of the 
field strength over the 2-dimensional boundary <9f2 3 of the 3-dimensional 
domain f2 3 . The analogous is true for the flux integral over dCl 3 . 

and, since the volume is arbitrary, we have the local version of magnetic flux 
conservation as 

d {i F jk] = 0. (31) 

We substitute the decomposition (|23|) into (|31|). Then we find the homo- 
geneous Maxwell equations, 

d a B a = , e abc d b E c + d a B a = (32) 

or, symbolically, 

div£ = 0, cmlE + B = 0. (33) 

Thus both, the sourcelessness of B a and the Faraday induction law follow 
from magnetic flux conservation. Both laws are experimentally very well 
verified and, in turn, strongly support the axiom of the conservation of the 
magnetic flux. 

The recognition that Maxwell's theory, besides on charge conservation, is 
based on magnetic flux conservation, sheds new light on the possible existence 
of magnetic monopoles. First of all, careful search for them has not lead to 
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any signature of their possible existence, see 0. Furthermore, magnetic flux 
conservation would be violated if we postulated the existence of a current 
on the right hand side of fl3T|). Now, Eq.(|9]) is the analog of (|3T[), at least 
in our axiomatic set-up. Why should we believe in charge conservation any 
longer if we gave up magnetic flux conservation? Accordingly, we assume - 
in contrast to most elementary particle physicists, see Cheng & Li |3|] - that 
in Maxwell's theory proper there is no place for a magnetic current]] on the 
right hand side of (|3TD. 



7 Constitutive law (axiom 4) 

The Maxwell equations ( |14|) and ([U]) or, in the decomposed version, (|T9|) and 
(P3|), respectively, encompass altogether 6 partial differential equations with 
a first order time derivative (the 2 remaining equations can be understood as 
constraints to the initial configuration). Since excitations and field strengths 
add up to 6 + 6 = 12 independent components, certainly the Maxwellian set is 
underdetermined with respect to the time propagation of the electromagnetic 
field. What we clearly need is a relation between the excitations and the 
field strengths. As we will see, these so-called constitutive equations require 
additional knowledge about the properties of spacetime whereas the Maxwell 
equations, as derived so far, are of universal validity as long as classical 
physics is a valid approximation. In particular, in the Riemannian space of 
Einstein's gravitational theory the Maxwell equations look just the same as 
in (|14|) and ([H]). There is no adaptation needed of any kind, see [11] . 



If we investigate macroscopic matter, one has to derive from the micro- 
scopic Maxwell equations by statistical procedures the macroscopic Maxwell 
equations. They are expected to have the same structure as the microscopic 
ones. But let us stay, for the time being, on the microscopic level. 

Then we can make an attempt with a linear constitutive relation between 
H ij and F kh 

H ij = 1 tfJU p M = 1 fx ijH Fh ) (34) 

with the tensor density x that is characteristic for the spacetime under 
consideration. We require [x] = 1, i.e., for the dimensionfull scalar factor fac- 
tor / we have [/] = q/cp — q 2 /h = C/(Vs) = A/V = 1/VL. The dimensionless 



4 This argument does not exclude that, for topological reasons, the integral in (29) could 
be non- vanishing, as in the case of a Dirac monopole with a string, see flq] . 
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"modulus" x^ kl i because of the antisymmetries of TC^ and F k[ , obeys 

x w = _^iki = _ x ijik (35) 

Moreover, if we assume the existence of a Lagrangian density for the electro- 
magnetic field £ ~ Ti^Fij, then we have additionally the symmetries 

x ijU = x kiij ? x [ijki] = Q _ ^ 

The vanishing of the totally antisymmetric part comes about since the cor- 
responding Euler-Lagrange derivative of £ with respect to the 4-potential 
Ai identically vanishes; here = 2<9[jA,]. For \ % i M \ this leaves 20 indepen- 
dent components^. One can take such moduli, if applied on a macrophysical 
scale, for describing the electromagnetic properties of anisotropic crystals, 
e.g.. Then also non-linear (for ferromagnetism) and spatially non-local con- 
stitutive laws are in use. 

The simplest linear law is expected to be valid in vacuum. Classically, the 
vacuum of spacetime is described by its metric tensor g 13 = g 3% that deter- 
mines the temporal and spatial distances of neighboring events. Considering 
the symmetry properties of the density x ? the only ansatz possible, up to 
an arbitrary constant, seems to be 



X»* = V-detSmn • (37) 



Note that x^ kl is invariant under a rescaling of the metric gij — > fl 2 gij, with 
an arbitrary function Q(x l ). Using this freedom, we can always normalize 
the determinant of the metric to 1. 

As an example, let us consider the flat spacetime metric of a Minkowski 
space in Minkowskian coordinates, 





c- 2 
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-1 / 



(38) 



If we substitute (^) into ( p7|) and, in turn, Eq.([J7]) and / = J e / fio into (p4|), 
then we eventually find the well-known vacuum ( "Lorentz aether" ) relations, 

n ij = </^VV Fhi or V = e E, H= (1/^)8. (39) 
V Mo 



5 With such a linear constitutive law it is even possible to derive, up to a conformal 
factor, a metric of spacetime, provided one makes one additional assumption, see M. 
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The law converts Maxwell's equations, for vacuum, into a system of 
differential equations with a well-determined initial value problem. 

Acknowledgments : We are grateful to Marc Toussaint for discussions on 
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A Four-dimensional calculus without metric 
and integrals 

In a 4-dimensional space, in which arbitrary coordinates x % are used, with 
2 = 0,1, 2, 3, one can define derivatives and integrals of suitable antisym- 
metric covariant tensors and antisymmetric contravariant tensor densities 
without the need of a metric. The tensors are used for representing inten- 
sive quantities (how strong?), the tensor densities for extensive (additive) 
quantities (how much?). The natural formalism for defining integrals in a 
coordinate invariant way is exterior calculus, see Frankel |J. However, we 
will use here tensor calculus, see Schouten |TJ]] and also Schrodinger |I5|, 



which is more widely known under physicists and engineers. 

Integration over 4-dimensional regions — scalar densities 

Consider a certain 4-dimensional region f^. Then a integral over Q4 is of 
the form 

/ Ad 4 S, (40) 
J sit 

where d A S := dx°dx l dx 2 dx 3 is the 4-volume element which is a scalar density 
of weight —1. We want this integral to be a scalar, i.e., that its value does 
not depend on the particular coordinates we use. Then the integrand A has 
to be a scalar density of weight +1. In other words, when using the tensor 
formalism, the natural quantity required to formulate an invariant integral 
over a 4-dimensional region is a scalar density of weight +1. 

Integration over 3-dimensional regions — vector densities 

Now we want to define invariant integrals over some 3-dimensional hyper- 
surface Q3 in a four- dimensional space which can be defined by the parame- 
terization x % = x t (y a ), a,b,c = 1, 2, 3, where y a are also arbitrary coordinates 
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on fi 3 . Then we call 

^■■■=^ a 4w a 4^ dy,dy2dy3 (4i) 

the 3-surface element on Q 3 . This quantity is constructed by using only 
objects that can be defined in a general 4-dimensional space without metric or 
connection. It can be constructed as soon as we specify the parameterization 
of Q 3 . Here e^i is the 4-dimensional Levi-Civita tensor density of weight —1 
and e abc the 3-dimensional Levi-Civita tensor density of weight +1 on Q 3 . 
Furthermore, this hypersurface element turns out to be a covector density of 
weight —1 with respect to 4-dimensional coordinate transformation. With 
this integration element to our disposal, the natural form of an invariant 
integral over f2 3 is 

/ A^Si. (42) 

Therefore, the natural object to be integrated over Q% in order to obtain an 
invariant result is a vector density of weight +1. 

Integration over 2-dimensional regions — covariant tensors or con- 
travariant tensor densities 

Analogously, we can parameterize a 2-dimensional region Q 2 by means of 
x % = x l (z a ), a,P = 1,2, where z a are arbitrary coordinates on Q 2 - Then we 
can immediately construct the following 2-surface element 

<es<r-= r m~ f ^d^\ (43) 

where e Q/3 is the Levi-Civita density of weight +1 on Q 2 - This surface element 
is an antisymmetric second order covariant tensor density of weight —1. Then 
an invariant integral is naturally defined as 

(44) 

with A %: > being an antisymmetric second order contravariant tensor density 
of weight +1. 

Alternatively, one can write the same integral in terms on an antisym- 
metric second order covariant tensor := ^eijkiA kl and an antisymmetric 
second order contravariant surface element 

^■^^-^e^dz'dz 2 , (45) 
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such that 



Jn 2 2 Jn 2 2 



(46) 



Since extensive quantities are represented by densities, we would take the 
first integral for them, whereas for intensive quantities the second integral 
should be used. Analogous considerations can be applied to (fit]) and (f42|). 

Stokes' theorem 

Stokes' theorem gives us as particular cases the following integral identi- 



ties (see [T3J p. 67 et seq.): 

/ {d l j l )d A s= f j l d 3 s t , 
I {d j n ij )d 3 s i = I \n^d 2 s l3 . 

Ja* Jan, 2 



(47) 
(48) 



B Decomposition of totally antisymmetric ten- 
sors into longitudinal and transversal pieces 

Here we provide the decomposition formulas for totally antisymmetric co- 
variant and contravariant tensors, which are the natural generalization of 
the decomposition of vectors and covectors. 

We start by considering an antisymmetric covariant tensor of rank p, 
namely t^...* . Its longitudinal and transversal components are given by 



U 



■ll...lp 



{p + l)L m [m U n ... lph (49) 



where ki := <9j<7, and L l ■ := n l kj. They fulfill the following properties: 



n n £ ix .., p = 



^ ^ii...i p ^ii...ip • 



(50) 



For p = 1, 2, 3, 4 we can explicitly write: 



p 


quantity 


definition 


explicitly 


1 




Li i u m 


ti kiU m 


2 








3 


±77 

l - y ijk 


3L m ^U rn \jk] 


n {kiUynjk -\- kjU m ki -\- k\JJ m %j) 


4 




AL m ^U m \ jk i] 


n {kiUmjki kjU m kn -\- kfoUfniij k]U m iyk) 
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Now we turn to U* 1 "^, an antisymmetric contravariant tensors of rank p. 
We define the decomposition as 

±yi 1 ...i p _ p ym\i 2 ...i p ] ^ yii-ip _ ^ _|_ l)/,l m yh---i P ] _ 

They fulfill 

kh yh-i P = o ? fcii ±yii...i P = ^yii-ip . ( 52 ) 

For p = 1, 2, 3, 4 we have the following explicit expressions for the longitudinal 
components: 



p 


quantity 


definition 


explicitly 


1 




U V m 


n l k m V m 


2 




21,1*1 


k m [n l V m3 - n 3 V mi ) 


3 


±yijk 


or[i| T/m|ifc] 
m v 


k m (n i V mik + n j V mki + n k V mij ) 


4 


Lyijkl 


AT[i\ ym\jkl] 


k m {n i V rnjkl - n j V mkli + n k V mlij - n l V mijk ^j 



An analogous scheme is valid for the corresponding densities. 
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